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Abstract 

Equations of motion corresponding to the Henon - Heiles system are considered. 
A method enabling one to find all elliptic solutions of an autonomous ordinary differ- 
ential equation or a system of autonomous ordinary differential equations is described. 
New families of elliptic solutions of a fourth-order equation related to the Henon - 
Heiles system are obtained. A classification of elliptic solutions up to the sixth order 
inclusively is presented. 

Henon - Heiles system, elliptic solutions, nonlinear ordinary differential equations 



1 Introduction 

The Henon - Heiles model belong to one of the most famous and frequently used models 
in astronomy and some other fields of physics. For example, the Henon - Heiles equations 
arise in description of a moving star in the axisymmetric potential of the galaxy. The Henon 
- Heiles model is defined by the following Hamiltonian [TJ 

H = X - (x 2 + y 2 + a x x 2 + a 2 y 2 ) + x 2 y - |y 3 . (1.1) 

Equations of motion corresponding to the Hamiltonian H are of the form 

x u = -Oiix - 2xy, 
Vtt = -at 2 y -x 2 + ay 2 

As is well known, in the following cases 

a = —1, «2 = aq; 

o = —6, 0(2 o.\ arbitrary; (1.3) 
a = —16, a2 = 16aq, 

the system (11.21) is integrable [2HE]. The values of the parameters given in (11.31) can be 
obtained by means of the Painieve methods. It is a remarkable fact that the integrable 
cases (11. 3p of the Henon - Heiles model are related to the stationary flows of the Sawada - 
Kotera equation, the fifth-order Korteweg - de Vries equation, and the Kaup - Kupershmidt 
equation [_5j. While scaling similarity solutions of these partial differential equations can be 
associated with time-dependent generalizations of the integrable Henon - Heiles Hamilto- 
nians p 7 ]. Backhand transformations for integrable cases are given in articles [HI HUH]. The 
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general solutions of the Henon - Heiles equations are found only in integrable cases [6j[9j[Tn] , 
while only several families of exact solutions are known in other cases [T3HT5] . A powerful 
algebro-geometric approach was used to find certain families of elliptic solutions arising in 
the Henon - Heiles model in integrable cases [H1[I2]- Some simply periodic and elliptic so- 
lutions were obtained in articles [T3HT5] . A systematic study of all the situations, when one 
can get exact solution of the Henon - Heiles system in explicit form, has not been undertaken 
yet. 

It is a distinguishing property of a wide class of algebraic autonomous nonlinear or- 
dinary differential equation equations that all their meromorphic solutions are in fact el- 
liptic, rational in cot bz (b is a constant) or rational in z with z being an independent 
variable [T6HT8] . Explicit expressions and a way one can construct them in closed form 
were presented in [TTHT9] . In this article we describe a method, which can be applied to 
construct all elliptic solutions arising in the Henon - Heiles model. We emphasize that this 
method enables one to find and classify elliptic solutions not only in integrable, but also in 
non-integrable cases. We use this approach to make a subsequent investigation of elliptic 
solutions of a fourth-order equation related to the Henon - Heiles system. Let us mention 
that despite the existence of the general solutions in the integrable cases it is not an easy 
task to find all reductions of known general solutions to elliptic functions. 

This article is organized as follows. In section [2] we give a detailed description of our 
approach and present a method enabling one to find any elliptic solutions of an autonomous 
ordinary differential equation. In section |3] we consider the fourth-order equation associated 
with the Henon - Heiles system and study the local behavior of its meromorphic solutions. In 
sections 4, 5, 6 we classify elliptic solutions of the equation in question restriction ourselves 
to the cases of elliptic solutions up to the sixth order inclusively. 

2 Method applied 

Let us take an algebraic autonomous nonlinear ordinary differential equation 

E[w(t)] = 0. (2.1) 

Our goal is to describe a method, which can be used to obtain all its elliptic solutions in 
explicit form. An elliptic function is defined as a meromorphic function periodic in two 
directions. Any elliptic function can be determined by its behavior in a parallelogram of 
periods. The number of poles in a parallelogram of periods, counting multiplicity, is called 
the order of an elliptic function. If w(t) is such an elliptic solution of equation (12. ip . then 
equation in question has the family of elliptic solutions w(t — to) with to being an arbitrary 
constant. Equation ( 12. ip necessary possesses an elliptic solution if it admits at least one 
Laurent expansion in a neighborhood of the pole t — to. Without loss of generality we shall 
build the Laurent series in a neighborhood of the point t — 

p oo 

™W = £lF + X^ fe > °<M< £ - (2-2) 

k=l k=0 

Here p is the order of the pole t = 0. The Painleve methods allow one to find all Laurent 
series satisfying equation (12. ip . The following proposition was given in article [T9] . 

Proposition. Suppose Laurent series (12. 2p with uniquely determined coefficients sat- 
isfies equation (12. ip . - then this equation admits at most one meromorphic solution having a 
pole t = with Laurent series (12. 2p . 
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Thus we see that equation ( 12. ip may have at most one elliptic solution possessing the pole 
t = with Laurent series (12. 2p provided that all coefficients in ( 12. 2 p are uniquely determined. 

Our algorithm for finding elliptic solutions of equation ( 12. ip in closed form is the following. 
Note that we omit arbitrary constant to- 

Step 1. Perform local singularity analysis around movable singular points for solutions 
of equation (12. ip . 

Step 2. Select the order M of an elliptic solution w(t) and take K distinct Laurent series 

Pi (i) oo 

^HE^ + E^ °<l*l< £ - i = l,...,K (2.3) 

k=l k=0 

from those, found at step 1, in such a way that the following conditions 

f>«=0, f> = ^- (2-4) 

i=l i=\ 

hold. 

Step 3. Define the general expression for the elliptic solution w(t) possessing K poles ai, 
. . ., ax in a parallelogram of periods in such a way that the Laurent series in a neighborhood 
of the point t = is w^'(t — a*), % — 1, . . ., K (see [IZ1[T8]). In other words, take the 
following expression for w(t) 

K ( K pi /_-.\* (») jk-2 ) 

w(t) = e c»c(* -«*)+ EE ik^wh^ p{t ~ ai) + ~ ho - (2 - 5) 

i=l I i=l fc=2 ^ J 

Here p(t) is the Weierstrass elliptic function satisfying the equation 

(p t ) 2 = 4p" - g 2 p - g 3 , (2.6) 

C(t) is the Weierstrass ^-function, h Q is a constant. 

Step 4- Find the Laurent series for w(t) given by (12. 5p around its poles a\, . . ., ax- 

def def 

Without loss of generality set a\ = 0. Introduces notation A t = p(a,i), Bi = pt{a>i), i — 2, 
. . ., K, and 

K 

h = h - E^iC(^) - E c -2^)' ( 2 - 7 ) 

Using addition formulae for the functions p and ( (see [T7l[T8] and equalities (12.91) below) 
rewrite expression (12.51) as 



»(*)= EE 



(A; - 1)! rft fc " 2 \ 4 

i=2 k=2 v 7 1 \ 



Pt(t) + BS ' 2 
. pit) ~ A 



.Step 5. Require the Laurent series found at the fourth step to coincide with the corre- 
sponding Laurent series taken at the second stem. Form a system of algebraic equations. 
Add to this system the equations Bf = 4Af — g 2 Ai — g 3 , i = 2, . . ., K. The number of 



3 



equations in the system should be slightly more than the number of parameters of elliptic 
solution (12. 8p and equation (12. ip . Solve the algebraic system for the parameters of the ellip- 
tic solution w(t), i.e. find h , g 2 , #3, A4, Bi, i = 2, . . ., K. In addition correlations for the 
parameters of equation (12. ip may arise. If this system is inconsistent, then equation (12. ip 
does not possess elliptic solutions with supposed Laurent expansions around poles. 

Step 6. Check-up obtained solutions, substituting them into the original equation. 

Further, we note that equation (12. ip may possess the Laurent series with arbitrary coef- 
ficients. If such series are taken at the second step, then it is convenient to add the arbitrary 
coefficients to the list of parameters. In addition in the fifth step of the method the algebraic 
system can be also obtained in the following way: one substitutes all the Laurent series found 
in the fourth step into the original equation and sets to zero the coefficients at negative and 
zero powers of (t — a,) in the resulting equality. This approach uses the fact that an elliptic 
function without poles is a constant. Suppose expression (12. 5p is built on the basis of the 
Laurent series involving those that contain arbitrary coefficients; then one should verify that 
all the series are in fact distinct, especially if one forms the algebraic system with the help 
of the aforementioned approach. 

In order to find the Laurent series at the fourth step of the method we use addition 
theorems for the functions £ and p whenever the elliptic solution w (t) possesses more then 
two distinct poles inside a parallelogram of periods. In our notation the formulae given by 
the addition theorems can be written as 



((di - a,j) = C(cti) - CKO + 



Bj + Bj 
2{Ai - A 
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(Bi + Bj, 



p{a, t - aj ) = -A - Aj + 4 V (i4 ._^ )2 (2.9) 



p t (di - = -Bi 



B i + B ] ){l2A 2 i -g 2 ) BiiBi + Bj 



" KA-Atf ZiA-AjY 

The values C( a i); i — 2, . . ., K disappear from the resulting series provided that we use ho 
instead of h in the series and expression (I2.8p . 

Further, let us obtain the addition formulae for the following situation jp(aj) = p(aj), 
a,i 7^ aj. The Weierstrass elliptic function p(t) takes each value twice in a parallelogram 
of periods. Suppose 2u%, 2co> 2 are the basic periods of the Weierstrass elliptic function. We 
denote by U3 the sum w\ + uj<i- The double points are t = 0, t = cui, I = 1, 2, 3 with t = 
being the double pole. Since the Weierstrass function is even, we see that the following 
relation is valid p(2ui — t) — p(t), I = 1, 2, 3. Consequently, in the case jo(oj) = we 
obtain pt{at) = —p t (aj) and p(a, — %•) = p(2aj). Note that the points a i: aj are simple, this 
yields pt{di) 7^ 0, pt(a>j) 7^ 0. Applying the L'Hopital's rule in (12.91) . we get the addition 
formulae for the case p{aA = 

C(2oi) = 2C(a,) + 



AB 1 
(12A?-g 2 ) ; 

3(12i4?-0O4 (12,4? -g 2 y 



p(2a t ) = -2A t+ ' ^ p2 y2 > (2.10) 

9/10 A2 

p t (2 ai ) = -Bi + 



2Bi 32Bf 

While solving the algebraic system we can simplify computations if we find the values £(cij 
aj), p{ai — aj), ptifli — ctj) first and only then use addition formulae 02. 9p or (I2.10p . 



4 



With the help of our method one can construct any elliptic solution of equation (12. ip . 
Note that if equation (12. ip possesses only N distinct Laurent series in a neighborhood of 
poles, then the orders of its elliptic solutions are not more than J2f =l Pi, where pi (i = 1, . . ., 
N) are the orders of poles given by the local singularity analysis [T7J[rS]. Thus we see that 
our approach may be used, if one needs to classify families of elliptic solutions satisfying 
equation (12.11) . Further, we mention that in the case g\ — 27 g 2 = the elliptic function p{t) 
degenerates and consequently elliptic solution (12.81) degenerates. 



3 Local singularity analysis for the Henon — Heiles 
model 

The method presented in section [2] can be easily extended to a system of autonomous 
nonlinear ordinary differential equations. Nevertheless in this article we shall consider the 
forth-order equation 

20 

y U tt-2{a-A)yy tt -2(a+l)y 2 t -—ay i +(Aa 1 +a2)ytt+{Qa2-^a l )y 2 +Aa 1 aiy+AH = (3.1) 

satisfied by the function y(t). An advantage of this approach lies in the fact that the only 
function supposed to be elliptic is y(t), while this may not be true for x(t). Equation (13. ip 
can be obtained in the following way. We differentiate the second relation in (II .2p twice 
and use expression (II .ip and the first relation in (ll.2p to eliminate x 2 and x u - Note that H 
in (13. ip is the energy of the system. Further, let us suppose that the variable t is complex 
and a ^ 0. 

We use the Painleve methods to obtain the Laurent series satisfying equation (13. ip . The 
dominant behaviors and the Fuchs indices are the following 

y (1) (t) = -4, J = "l, 10, | ±^-23 - 24(7, 

6 r-^s (3 - 2) 

Note that we omit the arbitrary constant to- Consequently, equation (13 .ip may admit two 
families of Laurent series: 



£i) 

^H^+E4^ 

k=0 

If the parameter a is such that the first family does not have positive integer Fuchs indices 
with the exception of j = 10, then the series y^\t) in fact exists and possesses one arbitrary 
coefficient (in addition to the arbitrary constant to)- A similar statement is true for the 
second family. The Laurent series y^{t) with arbitrary coefficient Cg exists provided that 
the parameter a is chosen in such a way that there are no positive integer Fuchs indices with 
the exception of j = 5. In other cases additional arbitrary coefficients may enter the Laurent 
series. Thus such cases should be considered separately. 
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We would like to note that the equalities c^ 1 ] = 0, = are valid unless the corre- 
sponding series posses the Fuchs index j — 1. For the series y^(t) this situation takes place 
in the case a = —4/3 and for the series y^ 2 \t) such a situation occurs if a = —16/5. 

Let us take a = —4/3, then the Laurent series y^'(t) possesses two arbitrary coefficients 
c 2 > c s corresponding to the Fuchs indices j = 4, j = 10 accordingly. However this series 
satisfies equation (13.1 p under the condition 

2 72 
cti + — (4«i - 3a 2 )cti + — — (2a x - a 2 )(ax - a 2 ) = 0, c_i = cr}. (3.4) 
it ly^o 

Analogously, setting a = —16/5, we obtain the Laurent series y^ 2 \t), which may have two 

(2) (2) 

arbitrary coefficients c 3 , c 7 associated with the Fuchs indices j = 5, j = 9 accordingly. 
This series satisfies equation (13. ip provided that one of the following conditions 

(i) c_i = 0, c 3 (lQa 1 - a 2 )(16tti - 9a 2 ) = 0, c_i = ci|, c 3 = 4 ; 

, , , 675 (16«i - 5a 2 ) , 10125(16«i - a 2 )(16ai - 9a 2 ) (2) ( 3 - 5 ) 

m) c 1 H c 1 H = 0, c_i = c 1. 

1 ; ~ x 47872 ~ x 343146496 ' 1 ~ x 

holds. From (13. 5 p we see that the coefficient c\ is always arbitrary and the coefficient c 3 
is arbitrary either in the case (ii) or in the case (ii) with a 2 = 16a\ or a 2 = 16«i/9. 

Further let us consider the values of a, which correspond to integrable cases (ll.3p . We 
begin with a = —1. The Laurent series y^(t) exists only in the case a 2 = «i and possesses 
three arbitrary coefficients cjp, , c§ related to the Fuchs indices j = 2, j = 3, j = 10 
accordingly. The Laurent series y^ 2 \t) satisfies equation (13.11) under the condition 

c 3 («i-a 2 ) = 0, c 3 = c {2) (3.6) 
(2) 

and possesses at least one arbitrary coefficient q . In the integrable case a 2 = a 2 the series 

y^ 2 \t) admits another arbitrary coefficient c 3 2 \ Interestingly that the series y^(t) may exist 

(2) 

in the case a 2 7^ a 2 , but with the loss of the arbitrary coefficient c 3 , which is zero. 

In the case a = — 6 the Laurent series y^(t), y^ 2 \t) satisfy equation (13. ip at any values 
of the parameters cti, a 2 . The series y^ 2 \t) possesses three arbitrary coefficients: cf\ 
Cg 2 \ while the series y^(t) only two: Cg , Cg . 

Finally in the case a = — 16 equation (13. ip admits the Laurent series y^(t) whenever 
one of the following conditions hold 

(i) a 2 = 16oti, 

, . .. [14a;ia; 2 — 112a? — a-lA m 

{U) C8 = 47567520 ' C§ " C§ ' ( 3 ' 7 ) 

7658 o 3829 2 1057 2 71 s TT 

A = «i «i«2 aia, H a 9 + ii 

22275 1 59400 1 118800 2 178200 2 

In the case (i) the coefficients c 8 , arbitrary, while in the case (ii) there exists only 

one arbitrary coefficient: c$ . The Laurent series y^ 2 \t) satisfies equation (13. ip under the 
conditions 

(i) a 2 = 16«i, 

(U) cS 2) = o, 4 2) = o (3 ' 8) 
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and possesses three arbitrary constants q , c 3 , c 5 in the case (i) and only one q in 
the case (ii) . Again we note that one can construct the Laurent series in the situation 
a 2 7^ I6a\. 

In addition let us note that in the case a = — 2 the series y^'(t), y^ 2 \t) coincide. Such a 
series possesses an arbitrary coefficient c§ and exists under the condition C3 = 0. 

The local analysis shows that elliptic solutions of equation (13 .ip possess even orders. 
In the next sections we shall classify second-order, fourth-order, and sixth-order elliptic 
solutions. 



4 Second— order elliptic solutions of equation ( 13.11) 

In this section let us construct second-order elliptic solutions of equation fl3.ll) . Using 
results of section [3] and relation (I2.5p . we find the general expression for second-order elliptic 
functions that satisfy equation (13.11) : 

y(t) = c- 2 p(t, g 2 , g 3 ) + h . (4.1) 

Again we omit the arbitrary constant to- Constructing the Laurent series for this function 
in a neighborhood of its pole t = 0, we obtain 

V(t) = C -f + h* + |c- 2 t 2 + |c_ 2 t 4 + . . . . (4.2) 

According to the method of section [2] we should either compare the coefficient of the 
series found by the methods of local analysis with coefficients of series (14 .2p or substitute 
series (14. 2 p into the original equation and set to zero the coefficients at negative and zero 
powers of t. Since series (I3.3P may possess arbitrary coefficients, we shall use the second 
alternative. 

Let us present our results for the case c_2 = —3. We have found five families of elliptic 
solutions with the local behavior y^\t). The parameters ho, g 2 , 93 for the first family of 
elliptic solutions are the following 



a 2 -(a + 2)a 1 (a + 2)(3a - 2)(aa l + 2a 2 )a l + (7a + \2)a\ 

"0 — 77 — , 92 



93 



4(a + l) ' ^ 12(^ + 1)2(3^ + 4) 

(a + 2) aia 2 6i + S 2 - ASH (3 a + 4) (a + l) 3 



(4.3) 



(4.4) 



216(a+ l) 3 (3(T + 4)((T + 3) 
where we have used the designations 

8 X = 3 (3a 3 + 7 a 2 + 28 a + 4) a x - 3(7a 2 + 17a - 10)a 2 , 
5 2 = (a + 2) (3a 3 + 7a 2 + 28 a + 4) aa\ - (3a + 8) (5a + 9) a\ 

The first family exists provided that a ^ —1, a —3, a ^ —4/3 and appeared in article 
The parameters g 2 , g% for the second family are given by 

_ 4 {5h 2 + 5ai/i + a 2 } 320/ig + 480aifc§ + 222a 2 /i + 30a? + 12H ^ 
92 3 , 93 - 108 1 j 

This family satisfies equation (13. ip under the conditions a — — 1, a 2 — cti and possess the 
arbitrary constant h . In the case a = —3 we obtain the third family with the parameters 

, cki + a 2 9aa — 22a 2 ai + 33a? a ^ 
ho = g— , 92 = ^ • (4-6) 
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This family exists under the condition 

_ (a 2 - «i)(Q2 - 7ai)(a 2 + 7oi) ^ 
320 

and possesses the arbitrary constant g 3 . In the case a = —4/3 we obtain two families (the 
fourth family and the fifth family in our numeration) of elliptic solutions. The fourth family 
satisfies equation (13. ip under the conditions a 2 = «i (see (13.41) with c\ = 0). The parameters 
ho, (73 are the following 

7 a i 31 o 7 2 , 

4 ' y3 2160 1 60 y 15 V ; 

The fifth family solve equation (13. ip provided that a 2 = 2a,\ (see (13. 4p with c\ = 0). The 
parameters /io, #3 are given by 

7 2.2 
fco = -ai, 9s = ^192 ~ T^a 3 - —H. (4.9) 
30 135 15 

For the last two families the parameter g 2 is arbitrary. 

Further we proceed to the case c_ 2 = 6/ a. We have obtained two families of elliptic 
solutions with the local behavior y^ 2 \t). The parameters ho, g 2 , g% for the first family are 

a 2 a| a| u iq\ 

~ 2a' ^ 2 ~ 12' ^ 3 ~ 216 ~ ^8~' 

If a = —6 we can construct two-parametric elliptic solution, which generalizes (I4.10p . Thus, 
the second family exists in the case a = —6 and possesses the arbitrary parameter h . The 
other parameters are given by 

g 2 = -Qh {l0h + 4«i + a 2 } - 1a x a 2 , g 3 = 14/iq{20/i + 3a 2 + 12«i} 

hn ceiUo (4-11) 

+ y {4«x + 3a 2 }{12 ai + a 2 } + -^{^ + a 2 } - 2H. 

The families of second-order elliptic solutions with a = —3, a = —3/4 seem to be new, 
while the families with a = —1, a = —6 are reductions of known general solutions in the 
integrable cases. Concluding this section we would like to note that elliptic solutions we 
have found degenerate if the following condition is valid: g\ — 21g\ = 0. 



5 Fourth— order elliptic solutions of equation (13.11) 



In this section our aim is to find fourth-order elliptic solutions of equation (13. ip . Using 
formula (12. 5p and the results of asymptotic analysis, we see that if such an elliptic solution 
exists, then it necessary possesses two double poles inside a parallelogram of periods. With- 
out loss of generality let us omit the arbitrary constant to an d suppose that these poles are 
t — 0, t — a. There are two possibilities. Equation (13. ip may admit fourth-order elliptic 
solutions possessing simultaneously poles with the Laurent series given by y^'(t), y( 2 \t — a). 
In addition there may exist fourth-order elliptic solutions possessing poles with the Laurent 
series of only one type. The latter situation may take place only if the series corresponding 
to the poles t — 0, t — a are in fact distinct. In the first situation we get the following 
expressions for the fourth-order elliptic solutions 

y(t) = c ( -lp(t, 92, 9s) + c^lpit - a, g 2 , g 3 ) + h , (5.1) 
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where h Q = h + c[_ 2 p(a, g 2 , g%). The Laurent series of this function in a neighborhood of the 
poles t — 0, t — a are of the form 

V(t) = ^ + ho + 2Ac^ - Bc%t + { (sA 2 - f) cS + |c«} t 2 + ... 

(2) 

*(*) = jp^ 2 + h + 2Ac^l + - a) + { - f) c« + §<£>} (t - a) 2 + . . . 

(5.2) 

def def 

Note that throughout this section we use notation A = p(a, #2, #3), -B = p t (a, g 2 , g^). In 
the second situation the general expression for the fourth-order elliptic solutions can be 
written as 

y(t) = c_ 2 p(t, #2, ga) + K(t, 02, 93) + c- 2 p(t - a, g 2 , g 3 ) - b((t - a, g 2 , g 3 ) + h , (5.3) 

where h = h + c_ 2 p(a, g 2 , g 3 ) — b((a, g 2 , g 3 ), c_ 2 = cl 2 or c -2 = c -l- Finding the Laurent 
series of the function (15.31) in a neighborhood of the poles t — 0, t — a, we obtain the series 

V(t) = °-f + - t + h o + 2^c_ 2 + {bA - Bc_ 2 }t + I (3A 2 - |) c {2) 2 - h -B^ t 2 + . . . 

y(t) = , °~\ 2 - 73— + /i + 2Ac_ 2 - {bA - Bc_ 2 }(t - a) (5.4) 

[~b (X J V Ob 

+ {(3A»- !)<=.„- (*-„)» + . .. 

We see that the series given in (15.41) are distinct in the following cases: b 7^ or b = 0, B 7^ 0. 
Thus we should consider the values of a giving the Fuchs indices j = 1 or (and) j = 3 for at 
least one of the series y^(t), y^ 2 \t) (see table 15 



Table 5.1: Values of a with the Fuchs indices j = 1 or (and) j = 3. 



Values of a 


Fuchs indices 


Type of series 


4 

3 


-1,1,4,10 


yW(t) 


-1 


-1,2,3,10 


yW(t) 


16 


-1,1,5,9 


v w (t) 


— 16 


-1,3,5,7 


y (2) (t) 



Let us present our results. In order to find elliptic solutions of the form (15.11) or (15. 3p we 
calculate five coefficients Ck, k — 0, . . ., 4 in each of the series (I5.2p or (15 .4p . Substituting 
these series into equation (13.11) . we get 10 nontrivial algebraic equations. Further we solve 
the algebraic system with an additional equation B 2 = 4A 3 — g 2 A — g 3 . There exists a 
number of solutions (15. ip corresponding to degenerate elliptic solutions. We do not present 
these solutions here. The only case, which leads to non-degenerate elliptic solutions of the 
form (I5.ip . is a = —6. This family of elliptic solutions can be written as 

-Mt\92,93)+2A-^-^-, (5.5) 



»W = "4 



pt{t; 92,93) 
p(t]92,93) - A. 
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Table 5.2: Coefficients of relations (I5.10p . ( 15. lip in the case a = 



-4/3. 



7l 

72 
73 



18 V + 282 



jf^l ~r -YfCt20ti — 2772^2 



3923 n 2 



9 I 891 1 ^ 1188 2 J ~ 294030 2 1 



4 rr _ 

135 27225 



22552^,3 i 239873^,3 



a? + 



653400 



a 



o _ 1391287 2 , uuoo n i. n 
2 980100 "2"! ^ 3025 1 2 



310253 2 _ 604523 2 
392040 2 441045 1 



5638 ^,2, 



/^2 — 24 ft 2 — TgCH — 24^ 

2581 {3a 2 -4 ai }6 2 2 



/'i 



14256 



_ f 10961063 
^° 1 82328400 2 1 



2921 



54450 1 



a? + 



2921 
24200 



51975 

10961063 2 _ 
123492600 1 439084800 2 



11507^,2 

7 UL n 



"l 17325 2 1 46200^2 

22367087 „,2~l 1,2 _ 1439023 2 , 247217 3 
J 15681600 2 a l "r 10454400 2 



a 2 «i - 5io 



where the invariants g 2 , g-s are given by 

a 2 3 o ^ ^2 

go = 30A a 9 H aocn Oi, = 

y2 280 2 35 35 l ' y3 280 



(7280A 2 - 3a 2 2 + 16a 2 ai - 48a 2 ) (5.6) 



and the parameter A is one of the roots of the following cubic equation 

■ :> x H a 3 a?a 2 a 2 a^ 2a\ 



a 2 a\ 
357 



« 2 



1904 ~ 119 / + 918 + 367200 



+ 



11475 45900 11475 



(5.7) 



Note that we have used an addition formula for the Weierstrass elliptic function (see ( 12. 8p ) 
in order to rewrite expression (15. ip in the form (15. 5p . 

Now let us find fourth-order elliptic solutions of the form (15. 3p . We begin with the case 
a = —4/3. From table I5TT1 it follows that we should take c_ 2 = —3 in expression (15.31) . The 
local singularity analysis (see ( 13. 4p ) shows that the parameter b satisfies the relation 



2 72 
b + — (4qi - 3ct 2 )6 2 + — — (2ai - a 2 )(a 1 - a 2 ) 

11 LvZo 



0. 



Solving 11 algebraic equations, we obtain the following family of elliptic solutions 



y(t) 



pt(t; 92, gs) + B 



4 L P(t;92,93) - A _ 
where the parameters h , B, g 2 , gz are given by 



Pt{t\ 92, 9z) + B 
p(t]92,93) - A _ 



+ h, 



(5.8) 



(5.9) 



1 29 ,2 
ho = 6A - -a 2 + -ax + —b , 



B = -— (12A - 346 2 + 15a 2 - 20ai) 



g 2 = 20A 2 + 
93 = 12b 2 - 



3605 



3605 



35 



1188 
A 



a 2 



-ai 





4 


b 2 + 5 


a 2 ai 




3 1 



A b 2 + 5 a 2 cti A + 71, 



1188 



891 3 

(924{4«i - 3a 2 }A + 3081a 2 . - 8088a 2 a! + 5392a 2 ) + 73 



(5.10) 
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The parameter A is one of the roots of the cubic equation 

A 3 + /j 2 A 2 + [i\A + /i = 0, (5.11) 

while the parameter b 7^ is a root of equation (15.81) . The coefficients 71, 72, fio, fix, \i 2 in 
expressions (15.101) . (15. lip are presented in table 15.21 

Further we proceed to the case a = —16/5. Using table 15.11 we set c_ 2 = —15/8 in 
expression (15.31) . While solving the algebraic system, we obtain that from b = it follows 
B = 0, Consequently, we suppose that b 7^ and consider the case (ii) in relation (13.51) with 
c_i = b. We find the following family of elliptic solutions 
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2 



p(t;92,93) - A 
where the parameters ho, B, g 2 , g% take the form 



Pt{t\ 9i, 9z) + # 
p{t]92,93) - A 



+ h , (5.12) 



62, 2 5 15 „ „ / 2 5056 l2 5 8 , 

h = — b a 2 H A, 5 = -a 2 H 6 a 2 A 

45 32 4 ' V9 3375 72 15 

/1697 27152 32 \ l2 ir/l2 . 

a 2 = «9 ati A I 6 + 15A + 5i, (5.13) 

y \5049 25245 1 45 / ' V ; 

c ,2 /^ 41 35 17915 2 827 2 \ . . 

03 = o 2 o + a 2 o.\ a 9 a, ] A + 

y V 8976 287232 2 1122 1 J s 

and the parameter A is a root of the cubic equation (15. lip with the coefficients /xq, /^i, A*2 
given in table 15.31 The parameter b 7^ is a root of the equation 

4 675 (16«i - 5a 2 ) , 2 10125(16ai — a 2 )(16o;i — 9a 2 ) _ (5 14) 

+ 47872 + 343146496 ~ 



and the coefficients /io, fii, ^2 are presented in table 15.31 

Now let us consider the case a = — 1. Setting c_ 2 = —3 in expression (I5.3p . we see that 
the series y^(t) exists only in the case a 2 = a±. Since b = 0, we have B 7^ 0. Solving the 
algebraic system, we find the following family of elliptic solutions 

3 J r 36p f (t;fl 2 ,g 3 ) + ^5184^3 - 108Aa 2 - 6af + T2H \ 2 , _ . 1 

~~ _ 7 S oa I U \ A\ f + 6A -7; a li {0.10) 

4 [ 36 [pit- g 2 , g 3 ) - A) J 2 

where the invariant g 2 , g% are given by 

q 2 = —a 2 , o 3 = —a? - — # (5.16) 
yi 12 yd 216 1 18 

and the parameter A is an arbitrary constant. 

Further it remains to study the case a = —16. We take c_ 2 = —3/8 in expression 

/r)\ 

(15. 3p and come to a conclusion that b = 0, B 7^ 0. It follows from B 7^ that q 7^ 0. 
Consequently, the series y^ 2 \t) exists only in the case a 2 = 16ai as given in (13. 8p . We find 
the following family of elliptic solutions 



3 9p t {t; g 2 , g 3 ) + yg6ff - 81^ + 108Aa 2 - 48a? ,31 
yU " 32 1 9(p{t;g 2 ,g 3 )-A) ( 4 2 ai ' 
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Table 5.3: Coefficients of relations in the case a = 



-16/5. 



X _ 3433 2 I 17165 _ 114673 2 

1 15708 1 ~ 125664 2 1 4021248 2 

r _ _896 / i2 , f 14693 _ 235088 1 <_ 80008 2 , 10001 _ 455047 2 

w 2 45^ \ 8415 42075 1 J 147015 1 29403 2 1 7840800 2 

X _ 307 3 I 921 2 I 159785 3 18113 2 44 ff 

3 2420 1 ^ 7744 1 2 ^ 53526528 2 619520 2 1 225 



,. _ 416^2 



ii = J- f 23450 , 75040 1 i2 , 6737 _ 165353 2 _ 6737 2 

"1 225 I 561 2 ' 561 "l/ (y ~ 125664" 2 " 1 20106240 2 78540 1 

_ 1 f 46282237 , 370257896 2 , 413300471 2 1 i2 , 24427 2 1441193 2 

^0 225 I 6609141 "2«1 T 33045705 "1 ~ r 352487520 2/ w ">" 2238016 1 2 537123840 1 2 

24427 3 1 4216603 3 4_ tt 

2098140 1 ' 15469166592 "2 225 



where the invariant g 2 , #3 take the form 

g 2 = l5A 2 -fal g 3 = 4Aa? - 10A 3 + - ~H (5.18) 

and the parameter A is an arbitrary constant. Note that we may take any sign at the square 
root in expressions f)5.15p . (15.171) . 

We have finished the classification of fourth-order elliptic solutions of equation (13.11) . The 
existence of fourth-order elliptic solutions with a = —4/3, a = —16/5 was claimed in [T5] . 
while in explicit form these solutions are presented for the first time. Elliptic solutions (I5.5p . 
(I5.15p . (I5.17P are reductions of hyperelliptic general solutions of equation (13. ip with a = —6, 
a = —1, cr = —16. The degeneracy condition for the elliptic solutions constructed in this 
section is g\ — 27 g\ = 0. 



6 Sixth— order elliptic solutions of equation ( 13.11 ) 



Let us classify sixth-order elliptic solutions of equation ( 13. ip . It follows from formula 
(12. 5p and the results of asymptotic analysis that such a solution possesses three double poles 
inside a parallelogram of periods. Without loss of generality we omit the arbitrary constant to 
and suppose that these poles are t — 0, t — ax, t — Below we use the following notation 
r\ = c_2, r 2 = c^- There are four possibilities. Equation (13.11) may admit sixth-order 
elliptic solutions possessing simultaneously poles with the Laurent series y( 2 \t), y^(t — a{), 
y^\t — 02)- These elliptic solutions are given by the relation 

y(t) = r 2 p{t) + np(t - a{) + bQ{t - ax) + rxp{t - a 2 ) - b((t - a 2 ) + h , (6.1) 

where h = h + rxp(ax) + r 2 p(a 2 ) + &C( a i) — K( a 2) and b = whenever a 7^ —4/3. 
Further equation (13.11) may have sixth-order elliptic solutions with the Laurent series y^\t), 
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y( 2 >(t — ai), y( 2 >(t — a 2 ) in a neighborhood the poles t = 0, t = a±, t = a 2 . We obtain the 
following expression 

y(t) = r x p{t) + r 2 p(t - a x ) + £»C(* - «i) + r 2p(£ - a 2 ) - b((t - a 2 ) + h , (6.2) 

where h = h + r 2 p(ai) + r 2 jp(a 2 ) + &C( a i) ~~ K( a 2) and 6 = whenever <x 7^ —16/5. In 
addition there may exist sixth-order elliptic solutions of equation ( 13.1ft with the Laurent 
series in a neighborhood of the poles of only the first or the second type. Such solutions 
have the expression 

y(t) = np(t) + &iC(t) + np(t - a x ) + b 2 ((t - a x ) + np(t - a 2 ) - (&i + b 2 )({t - a 2 ) 

(6.3) 

+h , h = h + rip(ai) + r;p(a 2 ) + 6 2 C(oi) - (61 + ^)C(«2), 

where i = 1 or i = 2. In the case i — 1 (i = 2) we see that 61 = and 62 = whenever 
tr ^ -4/3 (cr ^ -16/5). 

The Laurent series of the functions (16.1ft . (16. 2p in a neighborhood of the poles t = 0, 
t = ai, t = a 2 take the form 

?/(*) = p + 2r j {A 1 + A 2 ) +/i - {r,- (B, + B 2 ) + 6 (A, - A 2 )}t + . . . 

Vtt) = u ~ + + ^ + + { r * 5 i + r o P h2 + bT 1>2 } {t-ai) + ... ( 6 4) 

(t — a\) t — a± v ' 

y(0 = 77 ^ x 2 - 7— + ^1 + + {^2 - - Wi )2 } (t - a 2 ) + . . . , 

[t — a 2 ) t — a 2 

where V\ = rj(Ai + A 2 + I\ j2 ) — 6Ji j2 + /io- Along with this we take i = 2, j = 1 for the 
function ( 16. ip and and we set i = 1, j = 2 for the function ( I6.2p . Throughout this section 

we use notation J lj2 d = C(a\ - a 2 ) - ((ai) + C( a 2)> ^1,2 = - a 2 ), Pi ,2 = pt(ai - 02), 

Ai = p(ai), B\ = p t {ai), A 2 = p(a 2 ), B 2 = pt(a 2 ). The Laurent series for the function 
(16. 3p is the following 

V(t) = ^ + J + 2r i( A i + A 2) + h + {(b 1 + b 2 )A 2 - Mi - n (B 1 + B 2 )} t + . . . 

y{t) = , r * V2 + rj- + ^2 + {r^Bx + P lj2 ) + (61 + 6 2 )T 1)2 - Mi} (* - o x ) + . . . (6.5) 

y( f ) = 17~ — \2 - + ^ + {n(B 2 - P 1)2 ) - 6 2 T 1)2 - hA 2 } (t-a 2 ) + ..., 

(t - a 2 ) 1 t — a 2 

where v 2 = r i (2A l + A 2 + T lt2 )-(b 1 + b 2 )J lt2 + h and u 3 = r i {Ai + 2A- 2 +T Xi2 ) -& 2 Ji )2 + /i - In 
order to find sixth - order elliptic solutions explicitly we substitute the series (16. 4p or (16. 5p 
into the original equation and set to zero the coefficients at negative and zero powers of £ in 
the resulting expression. Note that £ = t for the first series in (16. 4p or ( 16. 5p . £ = t — a\ for 
the second series in (16.4ft or (16.51) . £ = t — a 2 for the third series in (16 ,4p or (16. 5p . Further 
one should check that all the Laurent series in a neighborhood of the poles t — 0, t — a±, 
t = a 2 are distinct. Consequently, in the case of elliptic solutions (16. 3p we should consider 
the possibilities given in table 16. ll Indeed, the series are distinct if b\ ^ or (and) b 2 7^ 0, 
otherwise supposing that the corresponding coefficients cq, c\ of the series (16. 5p with b\ = 0, 
b 2 = are equal, we see that all other coefficient at corresponding powers also coincide. 
Forming an algebraic system we need five coefficients in each of the series c^, k = 0, . . ., 
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4. As a result we obtain 15 nontrivial algebraic equations. In addition we include into 
the system the relations Bf = 4Af — g 2 A { — g 3 , i = 1, 2, P 2 2 = ATf 2 ~ 92^x,2 — 93 and 
expressions given by addition formulae (12. 9p or (12.101) . We use the latter expressions when 
most of other equations are solved. If we find that A x = A 2 then we take equalities (I2.10p . 
otherwise - (12. 9p . 

Table 6.1: Values of a with the Fuchs indices j = 1 or (and) j = 2 or (and) j = 3. 



Values of a 


Fuchs indices 


Type of series 


i 


4 

3 


-1, 1, 4, 10 


yW(t) 


1 


-1 


-1 ,2, 3, 10 


yW(t) 


1 


16 

s 


-1 ,1,5,9 


y (2) {t) 


2 


-16 


-1, 3, 5, 7 


v w (f) 


2 


-6 


-1,2, 5,8 


y (2) {t) 


2 



Some of obtained elliptic solutions are rather cumbersome that is why we present here 
not complete list of sixth-order elliptic solutions. In the case a = — 1, a 2 = ot\ there exist 
elliptic solutions of the form (16. ip . We set r\ = —3, r 2 = —6 and use addition formulae to 
obtain the following explicit expression 



y(t) = 6p(t;g 2 ,g 3 ) 



Pt{t\ 92, 9z) + B x 
p(t;g 2 ,g 3 ) - A 1 



Pt(t; g 2 , g 3 ) + B 2 
p(t]g 2 ,g 3 ) - A 2 



+ h, 



(6.6) 



where the parameter h , A x , A 2 , g 2 , g 3 , B\, and B 2 are given by 



h = 6A 2 



Oil 

T 



A 1 = 0, 



A, 



33«i 



a 1 



132 



H 



92 

B 2 = 



132' 
-B 1 



93 



H 

234 



O'i 



2808' 



(6.7) 



2808 234' 

Note that we may take any value of the square root in the expression for B\. Along with 
this we suppose that H 7^ if a\ = 0. 

Equation (13. ip possesses sixth-order elliptic solutions of the form (I6.2p . With the help 
of addition formulae we rewrite such a family of elliptic solutions as 



y(t) 



3(a + 4) 3 [{pfa g 2 , g 3 )} 2 + Bf] 

■p{t\ 92, 93) + — : 7T5 + ho- 



v [p(t; 92, 93) - A x ] 

Note that in this case A\ = A 2 . The parameter h , Ai, B 2 , g 2 , g 3 take the form 
tt-i 



(6.8) 



h 







{(a 2 + 24a + 32)k - a 4 - 27a 3 + 176a 2 + 384a + 512} 



A, 



a 1 



92 



{k + a 2 + 17a + 16} {a + 4}, B\ = {k + a 2 + 17a + 16} 

17a + 16}. 



(6.9) 



1 2 {a 2 + 32a+16}, 93 = ^^ + 



a 



Ylv\ 1 J ' y " 216^ 3 

This family of elliptic solutions satisfies equation ( 13. ip under two constraints on the param- 
eters of the original equation: 

a 1 " :; 



— {(a + 1) (a + 16) k + 181a 2 + 272a + 256} , H = {76800z/ 3 - 1891a 7 

V\ i.ZV-y ^6.10) 



a - 2 



-(181a 6 + 3830a 5 + 77440a 4 + 967960a 3 + 4542080a 2 + 6952576a + 3594240)/?} 
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In relations (I6.13p . ( I6.10P the parameters k, ui, z/ 2 , ^3 are given by 



K = ±Vl81a 2 + 272a + 256, v x = a 3 
u 2 = 2(a 2 + 17a + 16) k + a 4 + 34a 3 4 



121 



^3 



12800 



a 6 + 



6187 
640 



■a° + 



314131 
1920 



-a* + 



7115 



-a 



+ 34a 2 + 140a + 272, 

502a 2 + 816a + 512, 

365171 , 133556 

H <x H 

200 75 



-a + 



3744 



(6.11) 



In addition we suppose that V\ 7^ and a 7^ —2. If a = —2, then the corresponding 
coefficients of the Laurent series in a neighborhood of the points a±, a 2 coincide and solution 
(16. 6p is no longer a sixth-order elliptic function. 

Further let us take a = —6. We obtain elliptic solutions of the form ( 16 .3p with i — 2, 
r 2 = —1. Rewriting this function with the help of addition formulae ( I2.9p . we get 



y(t) = p(t;92,93) - 



Pt{t\ 92, 9z) + Bi 

p{t;92,g3) - Ai 



Pt{t\ 92, 9z) + B 2 

p{t;92,g3) - M 



(6.12) 



where the parameters ho, B 2 , g 2 , 93, B\ are given by 

4a 1 + a 2 



h = A 1 + A 2 - 



B\ = 4A 3 2 + 



93 : 

Aa>ia 2 
105 



20 
675 

gi 

140 35 



Bo 



92 



+ 



4a? 



aia 2 

675 
4a 2 



4a 2 a 2 
675 



+ 



Ao- 



Bi, 

' J 140 35 
2H 

5400 + 27 ' 

af aia 2 4a 2 a 2 



4a 2 ai 
105 



(6.13) 



675 675 



675 



_a|_ 
5400 



2iJ 
^7 



The parameter A 2 is arbitrary and the parameter A\ is the following 



A 1 = -±± 
2 



a 2 af 
560 + 35 



aia 2 3^2 



105 



4 



(6.14) 



Again we note that the sixth-order elliptic solutions obtained in this section degenerate 
if the following condition #| — 27g 2 = is valid. Family (I6.8P is given here for the first time, 
while families (16. 6p . (I6.12p are reductions of known general solutions of equation (13. ip with 
-1, a = -6. 



a 



7 Conclusion 

In this article we have studied the problem of finding exact elliptic solutions of a fourth- 
order ordinary differential equation arising in the Henon - Heiles model. We have given a 
detailed description of a method enabling one to find all the families of elliptic solutions of an 
autonomous ordinary differential equation or a system of autonomous ordinary differential 
equations. 

The fourth-order equation related the Henon - Heiles model admits elliptic solutions 
of even orders only. We have classified all the families of second-order, fourth-order, and 
sixth-order elliptic solutions of the equation in question. We have given explicit expressions 
for all the families of second-order and fourth-order elliptic solutions and for dome families 
of sixth-order elliptic solutions. 

The method described in this article generalizes several other methods, such as the Weier- 
strass elliptic-function method [2"3Tf2"5] . the Jacobi elliptic-function method [2"oTl2"8] and their 
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different extensions and modifications [29|l30] . Most of the methods do not alow one to ob- 
tain all elliptic solutions. The only exception is the method of Conte and Musette [2T||22]. 
Since there is no need to find and integrate any additional equation as in the method of 
Conte and Musette, it seems that the method used in this article is more simple in applica- 
tion. The main ideas of the method we have described in this article is the following. One 
uses the local singularity analysis in order to construct an explicit expression for a solution 
and then by direct calculations one finds all the parameters of the solution and constrains 
on the parameters of the original equation if any. Let us note that Vernov and Timoshkova 
considered the problem of finding fourth-order elliptic solutions of the equation in question 
in the cases o = —4/3, a = —16/5 [15]. They used the method developed by Conte and 
Musette and found a first - order equation satisfied by such families of elliptic solutions, but 
the did not present these solutions explicitly. 

In conclusion we would like to mention that our approach [T71CEE] may be used to find 
and to classify simply periodic solutions of the equation in question. This problem will be a 
topic of further investigations. 
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